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Abstract A new classical theory of gravitation within the framework of general relativity is presented. It is based 
on a matrix formulation of four-dimensional RiEMANN-spaces and uses no artificial fields or adjustable parameters. 
The geometrical stress-energy tensor is derived from a matrix-trace Lagrangian, which is not equivalent to the 
curvature scalar R. To enable a direct comparison with the EiNSTElN-theory a tetrad formalism is utilized, 
which shows similarities to teleparallel gravitation theories, but uses complex tetrads. Matrix theory might solve 
a 27- year-old, fundamental problem of those theories (Sect. 14. For the standard test cases (PPN scheme, 
SGHWARZSCHlLD-solution) no differences to the EiNSTElN-theory are found. However, the matrix theory exhibits 
novel, interesting vacuum solutions. 
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1 Introduction 

In the last decades a variety of new classical gravitation theories as alternatives to the EiNSTElN-theory were 
proposed [33 . This increased interest is particularly motivated by some new phenomena, which can only be 
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explained with some additional presumptions (e.g. galaxy rotation problem, Pioneer anomaly, accelerated Universe 
expansion). On the other hand, new and enhanced experimental possibilities allow, to test their predictions [421 
I12II11] with unthought precision. We want to mention here only as representatives the Brans-Dicke theory [3], 
as famous example of a scalar-tensor theory and MOND 3 , which is supposed to give an alternative to "dark 
matter" . A recent discussion of this can be found in 5 • 

In this paper a new general relativistic gravitation theory, titled "matrix theory", is presented. It is derived 
from a matrix-trace Lagrangian, similar to the well-known Einstein-Hilbert action, but based on matrix 
formulation of the four-dimensional Riemannian spacetime. 

Like Einstein's original theory (without the "cosmological constant") it contains no free, "adjustable" parameters, 
except the Newtonian constant of gravitation G. Also, it does not introduce new, artificial fields, like Brans- 
Dickes scalar-field or others in vector-tensor theories. 

To compare it with the EiNSTElN-theory of gravitation, we generalize this Lagrangian with tetrad formalism, 
so that it contains four real, constant parameters {a,b,c,d). Each parameter set then characterizes a different 
gravitation theory, and it is shown that also the EiNSTElN-theory belongs to this class of theories, esp. it is 
described by the parameters (a, b, c, d) = (1, — ^, —^,0), while the matrix-theory is defined with (1, —1, 0, — ^). 

Matrix theory uses complex tetrads, because general base matrices t^i can only be represented with such 
tetrads. This might look unfamiliar to some readers, but we consider this similar to the situation in quantum 
mechanics. There we have a complex (non-measurable) wave function and real observables. Here, the tetrads 
themselves are also not measurable, only the - by definition - real metric is measurable. Moreover, it shows, that 
all test cases computed here (sec.|4)|, which represent macroscopic matter (real, symmetric stress-energy tensor), 
have solutions with real tetrads (for the PPN-test in section 14.31 this holds up to the requested approximation 
order) . 

This tetrad formalism shows, that the matrix theory can be regarded as generalization of the " teleparallel" 
approach (also called "distant parallelism" or "absolute parallelism") of tetrad gravity. This is based on an idea of 
Einstein, which uses a non-symmetric "Weitzenbock" connection with vanishing curvature tensor but nonvanishing 
torsion, which is extensively discussed until today 36,39,341 1251115) . A comprehensive overview can be found in 
[31j and also [IS], where the gauge aspects of the theory are stressed. If we would consider only real tetrads, 
the resulting theory ( "RMT" , see sec. (6)1 would belong to the one-parameter class of teleparallel theories, which 
are experimentally viable [32ll31j . like the teleparallel equivalent of Einstein's GR (TEGR) [27]. The teleparallel 
theory also allows an alternative coframe representation, which is used in [171118) to derive a conserved energy- 
momentum current, completely similar to the Maxwell- Yang- Mills theory. 

However, the usage of complex tetrads, which are necessary to map arbitrary matrices, and the new matrix-trace 
Lagrangian, containing a parity violating term, exclude typical "unphysical" tetrad vacuum solutions, which 
prevent a profound interpretation of previous teleparallel theories. 

In this paper, we use the conventional Levi-Civita- (or Christoffel-) connection, which is built from the 
metric tensor (see eq. ()116p ff.). The tetrad formalism here serves only as a general mathematical tool to compare 
different theories. Instead of the tetrads, we consider the base matrices as the fundamental entities. Matter 
influence to geometry (field equations) is mediated via these r^^, resp. Pfj,u, but matter reacts to geometry only due 
to metric (equations of motion). Consequently, the geometrical stress-energy tensor is potentially not symmetric 
and real, but it is forced to be so, since it is equal to matter tensor (this is sometimes discussed differently for the 
teleparallel theory, see e.g. [36], p. 15). 

However, many of the general tetrad computations presented here, are mostly standard (esp. the representation 
of the RiCCl-scalar by tetrads and the derived stress-energy tensor) and can be found at various places and in 
various contexts. E.g. our equation ()55p is equivalent to eq. (1) in [25]. The reason to sketch them here nevertheless, 
is to give a homogenous presentation with consistent notations. This allows the reader to follow them without the 
necessity to check several sources with different names for the same variables. 

The readers will surely notice, that the hermitian matrices introduced in eq. ^ can also be regarded as second 
order WEYL-spinors (see e.g. [23], p. 59 ff) with respect to their LORENTZ-transformation rule. However, we do 
not discuss quantum mechanical effects or quantum fields here, to limit the extent of the paper. 



2 Formulation of the matrix theory 

To give a clearer picture, we start in l2.1l with the matrix representation of the widely known tetrad- (or " Vierbein"- 
) formalism of general relativity (e.g. [2H], [33)! which is given by hermitian matrices (real tetrads), and generalize 
this in section \2l2\ to general complex matrices (complex tetrads). 
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2.1 Hermitian matrix representation of tetrad formalism in relativity 

Here we want shortly sketch, how the main tensors and equations of general and special relativity can be 
represented with hermitian matrices. This representation does not offer new equations, but it needs less 
independent prerequisits (metric signature. Maxwell eq.), than the usual component formulation. As far as we 
know, this cannot be found in the literature in this compact form. 

Tetrads are four real, covariant spacetime vectors, which are defined in each point of the spacetime. We 
denote them here by e'^{x^), where a = 0, 1, 2, 3 is the tetrad index and = 0, 1, 2, 3 the spacetime index (in this 
paper Greek letters /i, v,a,l3, . . . are used for spacetime and Latin letters a, b,c,d, . . . for tetrad indices). One of 
the first physicists, who used them for GR (titled as "four- legs") was M0ller, see e.g. his basic paper 27,. As many 
others, he regards them as the fundamental gravitational field variables, instead of the metric g^u- Moreover, 
tetrads are also a useful tool in geodesic applications of general relativistic problems 22,. 

Each individual tetrad (denoted by a certain fixed "a") is a covariant tensor of first rank. When r^^f, = 
diag[l, —1, —1, —1] denotes the Minkowski- metric, the metric tensor g^i/ is expressed as 

a b 

Ot^iy = Vabef^eu- (1) 
Tetrad indices a,b, . . . can be shifted with rjab and rj""^ while spacetime indices fi^v, . . . axe shifted with g^i, and 

The contravariant (inverse) tetrads then fulfil two orthogonality relations, which are {5 is the usual 
KRONECKER-symbol) : 

6^2 = 5^: and e^^ = <5t (2) 

By contracting with or eJJ any spacetime index of any symbol (tensor or non-covariant entity), can be trans- 
formed into a tetrad index, and vice versa, e.g. 

A^'el = A'' ^ A^'^A^'e't. (3) 

Matrix representation: with the tetrads one can construct four complex, hermitian 2 x 2-matrices, using the 
generalized Pauli spin matrices uq = (J^) = I2, <y\ = (10)5 ^2 ~ (i 0)' ''"3 ~ (o-l)' which we will denote with 

Tfi = e^gg . (4) 

This definition is very similar to the expression of spinor components of tensors with the help of Infeld - van der 
Waerden symbols ( 33 , p. 123 and [33, p. 48) g^^' = -^^a^' where A,B' G [1,2] are the spinor indices. 

These four matrices are hermitian by construction, linearly independent, and can replace the tetrads, since 
eq. ((4]) is an invertible map. They form a basis in the vector space of 2 x 2 matrices, like the four crfe0 We will 
denote them here as "base matrices". Some general relations with these matrices are listed in the appendix lA.il 

Any tensor of first rank with contravariant components A^ can be expressed as hermitian matrix A by 
(boldface Latin letters A, B, . . . as well as Greek letters T,a,p, . . . shall denote hermitian 2 x 2-matrices here) 

A = A^T^. (5) 

Since the base matrices build a covariant "tensor- matrix" , this means that the matrix A is actually invariant 
under all transformations x'^[x'^). Of course, also the infinitesimal line element (1-form) can be expressed as the 
matrix dx = dx^Tfj, and transformation equations are 

dx'^ — ^—^dx'^ — a'^dx" and Tv = a^Jr^. (6) 

This transformation rule for the base matrices states, that all components are transformed with the same 
coefficients (like the tetrads). 

A novelty of the matrix notation, in contrast to usual tetrad notation, is, that it defines an inner product 
(matrix product) and with the help of this, the need to postulate a Minkowski norm for the tetrads (with its 
sign-arbitrariness) disappears. This hermitian matrix-algebra can be seen as special representation of Heestenes' 
"space-time algebra" (STA), which is widely discussed in the literature, especially for the DiRAC-theory, see e.g. 



^ More generally, any set of 4 hermitian matrices cr^ can be used as basis, that preserves the orthogonality \T{ct'^ct'„) = 
Vmn- This is in close relation to the transformations described in eq. II15I I. 
^ I.e. every (hermitian) matrix A can be expressed as linear combination A = a'^r^, with complex (real) coefficients a^*. 
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The norm of a tensor is the simple matrix-determinant 

|A| =.g^^A^yl'^ = A^^^. (7) 

This is easy to derive from the properties of the (7^, namely ^'T{ami7n) = Vmn, where T(A) denotes the trace and 
A the "adjuncted" matri>|f| of a matrix A. This simple norm definition is only possible for a four-dimensional 
RiEMANNlAN spacetime with Minkowskian signature — , — , — ]. 

Additionally we have to introduce the contravariant basis = g^'^Tu and can derive the orthogonality 
relations 

-T(rpr,.) = Qf,!, and -Tirf.f") = 5';^. (8) 

If the matrix theory is formulated without tetrads, the first equation is to interpret as the definition of the metric 
tensor g^i/ and the second as the definition of the matrices r^. 

The more general scalar product of two tensors A, B has a similar matrix representation like the norm in eq. 



ir(AB) = g^.A'^B'' = A^B^. (9) 
The inverse relation of eq. Q is the trace expression (always real) 

A^ = i7-(Ar^) and A^^^T{Af^). (10) 

Tensors of higher rank are expressed by sets of hermitian matrices, e.g. a general tensor of second rank with four 
matrices 

Ap = ^^^r'' ^ Af,^ = -T{Af,fu). (11) 

With the above definitions the complete apparatus of special and general relativity can be drawn in matrix form. 
E.g. the covariant derivative of the basis is computed like for a conventional vector 

The Christoffel symbols rj^i, defining the connection here, have to be derived metric compatible from eq. ([T| 
(see appendix [as]). The matrix-representation of the antisymmetric second covariant derivatives of the basis then 
gives a definition of a Riemann tensor matrix, which is very similar to the standard formula: 

Interested readers can also have a look at [19], where representations of main topics of special relativity (e.g. 
electromagnetism, DiRAC-equation) with matrices are shown. As single example we cite here the matrix repre- 
sentation of Maxwells equations {d cr^-^jr is the partial derivation operator matrix, F = {E^. + iBi^)a^ the 
trace-free electromagnetic field matrix (non-herm.) and J is the hermit, current matrix)lf| 

dY = 3. (14) 

This is only one matrix eq., but it contains 8 real (4 complex) component equations, which are the four 
homogeneous (anti-hermit, part) and the four inhomogeneous Maxwell eqs. (hermit, part). 

Additionally to local spacetime covariance, the matrix equations exhibit another, independent global symmetry: 
If all matrices are synchronously transformed with one constant, unimodular matrix T (i.e. jT| = 1), preserving 
their hermitian property: 

A TAT^ (15) 

then obviously all relations, e.g. the metric in eq. ([SJ, remain unchanged. The transformation matrix T then 
contains 6 real parameters and it is easy to show, that it can be identified with a LORENTZ-transformation in a 
local MiNKOWSKl-coordinate systemlf] Consequently, in those coordinate systems (locally) both transformations 
may be combined arbitrarily. 



This term is not widely used in English mathematical textbooks. We define it here as A = |A|A~^. Please note, 
that only for 2x2 matrices we have the linear map A A and consequently only then |A| = i7~(AA) is a bilinear form. 
See appendix I A. II for a more detailed discussion. 
* In flat Minkowski spacetime we use = cr^ = const. 

^ The group of matrices T with complex elements, satisfying |T| = 1, is commonly denoted as SL{2,C). It is the "double 
cover" of the Lorentz-group, because both matrices T and (— T) perform the same Minkowski-space rotation. 



Matrix theory of gravitation 



W. Kohlcr 5 



To describe curvature in Riemannian geometry, we define the "rho"-tensor-matrix, as the antisymmetric 
partial derivative of the basis 

PfJ,l^ — Tfj,.jy — T,y^^. (lb) 

The tensor property (covariant transformation rule) of this matrix-tensor is evident. It consists of 6 hermitian 
matrices and thus contains 4 x 6 = 24 real components. From /9^,y = follows the vanishing of the RiEMANN-tensor 
~ i^-S- by the derivations in the appendix IA.3|I i.e. the spacetime is flat. On the other hand, for a flat 
spacetime we can always find a coordinate system with = const, and consequently p^i, = 0. Due to the tensor 
property this equation remains true, if an arbitrary coordinate transformation is applied. 

From the definition of p^^ and the basis in eq. (|4]) we get the tetrad formula (we use the common []-bracket- 
notation, but omit a frequently used factor 1/2) 

Pa-r = (eS,7 - e^,a)'^x = e^^ ,^^ax, (17) 

where ej^ is the "nonholonomity" 15 . Since the Pauli matrices are constant, it is evident that p^jy is the 
matrix representation of the exterior derivative of the basis 1-forms 0^ — e^dx^ . 

The expressions are further simplified by transforming the spacetime indices a, 7 into tetrad-indices. For this 
purpose, we define the antisymmetric tetrad expressions r^c ~ —Tea by (Schouten, [3S], pp. 99, denotes them as 
"objects of anholonomy" f2ac)- 

rac = ej^^^je^ec'. (18) 

These 24 coefficients Vac can be classified into two types. For 12 of them the upper index x is equal to one of the 
lower. They will be denoted here as "r-doublets" . The other 12, where all three indices x ^ a ^ c are different, are 
denoted as "r-triplets". This classification is independent of the coordinate system, because the Vac are invariant 
under all coordinate transformations. 

By Cartan's first structural equations one can see, that these terms are closely related to the "Ricci rotation 
coefficients", which may be defined from the covariant derivative (also see appendix IA.3[) 

y-is ^fi/ u u s i/s.s6/ c. c w /1n^ 

Qmn = e'mene^-v = -{rmn + -n (llmcrnb + Vncrjnh))- (19) 

From these one can directly derive the tetrad representation of the curvature tensor (see eq. (|114p ff.) 

R^mnl = ep(<5fC/m„ - 5n5mi),A + GxyiSmT^i + SnGml " ^^Qmn)- (20) 



2.2 General matrices and complex tetrads 

For the matrix representation presented above, it looks straightforward, to consider general instead of special 
(hermitian) matrices as basis r^. Another motivation comes from quantum mechanics, which cannot be formu- 
lated without complex wave functions. Therefore one may hope, that the ideas presented here can help to find a 
new link between quantum mechanics and gravity. However, this is not the topic of this paper, which covers only 
classical gravity. 

On the other hand, for tetrad gravity in usual formulation, it makes no sense to introduce complex - instead of 
real - tetrads, because the field equations are not altered. This is only the case, if we use the matrix-LAGRANGlAN 
defined in 12.31 which has additional complex terms. 

For general base matrices, we have to generalize the metric definition in cq. ((Sjl, because the distance 
ds^ — g^udx^ dx'^ must always be a real quantity. 

Regarding that under coordinate transformations the hermitian conjugated matrices obey the same transfor- 
mation rule as in eq. Q (the transformation coefficients are real), the appropriate definition i^ 

3m- =^ -Jirlfu + rlf^,) = iKr(r^r.), (21) 

which is symmetric and real for arbitrary matrices and for hermitians r^^ = it is equal to the definition in 
It formally resembles definitions of quantum mechanical observables, e.g. the DiRAC current and is invariant 
under unitary U(l) (phase) transformations — >■ e^'^r^^, additionally to its T-invariance described in (|15|) . 

^ One might also discuss to use the complex value (without 5R), which would define a hermitian metric tensor = Qvfi. 
For all equations, where only the symmetric part of occurs, e.g. the equation of motion, it is equivalent. This form also 
allows the usual index shifting with g. 
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Also the scalar product of two tensor matrices A = a^Tfj,, B = b^Tfj. (with a^,b^ = real) is to define 
consistently as the real number 

(A • B) '^^^ iKr(A^B) = a^b'^^^Tirlf,) = a^fo^'p^.. (22) 

With this definition all equations of general relativity stay valid, except the field equations. For clarity we add, 
that for all physical problems covered here, we consider only strictly real and symmetric stress-energy tensors T'^'^ 
of matter. This requires, that all possible imaginary and anti-symmetric parts of the geometric tensor also vanish. 
We do not discuss possible implications of those terms, instead we demand that all must be zero, for all classical 
gravity problems in this paper. Then, e.g. for the test cases solved with real tetrads in section |4l the imaginary 
parts form additional constraints, compared to a corresponding real-tetrad theory ("RMT", see sec.|S}. 

For the definition of the contravariant base matrices r'' we cannot use the metric g'^'^ here anymore, but 
the second eq. of ((Sjl gives an unique definition. The contravariant transformation rule stays valid, due to this 
orthogonality relations. 

If we want to utilize the tetrad formalism for general bases, we have to use complex tetrads in the decomposition 
— e^^aa- The inverse tetrads are also to define with their orthogonality relations in eq. ([2}, index shifting 
with g^'^ is also not applicable for them. 

We have to add here, that for general (non-hermitian) matrices r^, the metric is not necessarily locally 
Lorentzian. However, this is always true for the physically important case, when the imaginary parts of the 
tetrads are small (e.g. for the PPN-tests in sec. I4.3|l . General matrices can be decomposed into a hermitian and 
anti-hermitian part, and correspondingly the complex tetrads ej^ into real and imaginary parts: e'^ — + ih'^ 
(with h^t = real). The metric definition (f2T|l then gives g^i, = (//f/^ 4- hf^h^)riab- 
It can be shown, that it is locally Lorentzian, if all imaginary parts are small: jj/i^jj <^ 1, Va,/i. 



2.3 Lagrangian of matrix-theory 

It is an important feature of the EiNSTElN-equation in general relativity, that it can be derived from a LAGRANGIAN 
£. (see e.g. [33), namely its geometrical part equals the curvature scalar Ce — R- 

For deriving the stress-energy tensor and the field equations one has to find the stationary solution of the action 
integral 

d'^x^/^Cigf^i,, g^ij^x) (23) 

by variation of the metric tensor Sg^m- The same holds for the matrix theory, where we postulate another scalar 
based on the "rho" -tensor- matrix defined in eq. (|16|) . As explained, this tensor-matrix also characterizes the 
curvature of spacetime and it is straight forward to construct a theory of gravity based on this tensor-matrix. 

Here we construct the "matrix-LAGRANGlAN" Cz as real, scalar, bilinear form from the matrices p^^v and p^^jy. 
We request the same symmetry as for the metric definition ([21}, i.e. global T-invariance forces, that the matrix 
factors in the trace must build a bar-alternating sequence. 

There exist only two distinct tensor matrices, that can be built by bar-alternating contraction of pfiu, namely 

def def 

Pfi — f'^ Pfiw and Vfi — PfiwT'^ ■ With the request of unitary (7(1) invariance we postulate the following La- 
GRANGIAN, which is also quadratic in the first derivatives^ 

(24) 



Cz \^{r-^fPpivp) 



This expression is a real function of the t^,!"^ and their first derivatives 'Cz(r^, r^, r^,iy, r^_^) and contains no 
adjustable parameters. By construction, it is invariant under arbitrary coordinate transformations and constant 
(global) T-transformations described in eq. (llSp . Considering its additional unitary invariance under — >■ e^'^r^, 
we find that the symmetry group is SL{2,C) x U(l), which is a supergroup of 517(2) x (7(1), the important group 
of standard electro- weak "GSW-theory" (see we only discuss global symmetry here). 

For completeness we have to add, that H24|) is of course not the only possible form. In general, every bar- 



alternating permutation of the 6 factors t'^\t'''\ p\^,T^ , ppg exhibits the same symmetries and its tetrad 
Lagrangian has the common form H32p . But if we request, that the contracted forms p^ and i;^ should occur, 
but no doubly contracted matrices (hke f'^p^uT^), then only four alternatives remain: KT(r"'ff'^Xc,^), where 



The more general LAGRANGIAN of "viable" theories Cv, which is discussed in scction[3]4]for comparison, can be written 
in the same form, with an extra term (exhibiting the same symmetries, but real by definition): Cv(c) = Cz + ^'T{t'^'^ Vj^ply), 
where "c" is a free, real constant. With the terminology of section [3] for the extra term holds \'T(t'^^ v pply.) = Cc — 2£i,. 
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^ap ~ Pa'^^li^ Pa'^p^ Pa^Jp, Pa^p, respectively. The fourth alternative gives a completely similar Lagrangian 
as the first eq. ()24|) (namely £ = C{1, —1,0, +^) in eq. ()32|) . i.e. only the odd parity term £i has opposite sign, 
which does not affect any conclusions), while the second and third form have no odd parity term. 

To derive the field eqs., similarly to above eq. ((23)) . one could vary the base matrices Sr^ instead of &?/jJfl 

I ^ J d^x\\r\\C{r^„T^,,^) ^ SI ^ J d^x\\r\\T{6r^,T'') (25) 

This derivation of the stress-energy tensor matrix would be straightforward. But instead of this, we give an 
equivalent derivation with the use of tetrads in the next paragraph. This has the advantage to be more general 
and so allows a direct comparison to the EiNSTElN-theory. 

Tetrad representation of Cz'- 

With the terms in eq. H17p the Lagrangian in (|24[) can be rewritten as (all rj^^ are scalar and can be drawn out 
of the matrix trace) : 



j?ltT{r'^^f^pUf''^Pf,sf') = j?lt{{r^acrrlT{a^a''a..a^aya^)) (26) 



The trace of 6 PAULl-matrices above is computed using the techniques in the appendix lA.ll If we define for 
abbreviation the two contracted terms {ra is constructed on only of r-doublets and t°' only of r-triplets): 

<ie/ X J .a d,ef 1 y .abed /r>-s 

Ta = rax and t = -riytr^^A , (27) 

where A"'^'^'^ is the completely antisymmetric symbol, with A^^"^^ — 1, the result is 

Cz = ??""(rmr* - r°„j,(r^a)*) +i (t°r* - f°*ra) . (28) 

Both terms Cr and iCi are evidently real (£* = —Ci). 

The explicit appearance of the imaginary unit "i" in this formula is a consequence of utilizing PAULl-matrices CTq 
as basis in = eJJ era If| 

We note, that Ci has "odd parity" (due to the factor A) in contrast to all other terms of Cr and eq. (I32p ) with 
respect to the tetrad space (the tetrad parity operation is equivalent to the matrix transformation — > f^, 
which inverts the three spatial tetrads e^,k = 1, 2, 3). 

We now discuss the implications of using real or complex tetrads for the variation principle. The variation 
gives the definition of the stress-energy tensor components by 

J^^-= J + iSC, =^ J Se'^T^ + {Se'^riT-^r = 0. (29) 

If we consider only a priori real tetrads, like in conventional tetrad theories, also the variations must be real, i.e. 
{Se^)* = Se^, and the variation principle gives only (T^)* + T'^ ~ 25R(r^) — 0, which means that Ci does not 
contribute in this case (the resulting theory "RMT" is discussed in secfS)). 

For the case of potentially complex tetrads, both variations Se^, (Se^)* are independent and we get the full complex 
eq. = and both Cr, iCi contribute to T^. Of course, then it is sufficient to consider only e.g. the variation 
of Se!i, because the second leads to the same eqs. 



* Here is to replace g = ||t||, where ||r|| is defined as absolute value of tlie determinant |t| of all 4 X 4 components 
of the basis. For the variation of this term one has to use 6|t| = ■!^T(5r^r'') 

^ They are well suited as basis for hermitian matrices, but not the best choice for arbitrary complex matrices. Another 
choice are the matrix-components themselves, which leads to a spinor-like notation = tjJ^i?AB, A,B = (1,2), with 

'^12 = (qJ) , . . .. The Lagrangian, computed with these 16 complex terms t^^, instead of the tetrads ej^, 
is somewhat simpler. However, they must be transformed into tetrads anyway, to describe local MiNKOWSKl-systems and 
for the test cases of section |4] 
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3 Generalized Lagrangian in tetrad-form 

To be able to make comparisons between all possible tetrad theories and to find a general expression of T^, we 
generalize the Lagrangian of eq. (|28p to a more general, real bilinear form of the r^^, with constant factors 



bed 



C = rl,{ryj* Hfy^'. (30) 

The expression C is real = for arbitray complex rj^^, if and only if H^^"'^ = {H^'t^ f holds. 

We discuss here a general, LORENTZ-invariant, bilinear forrrF^I that contains four free, constant parameters 

a, b, c, 43 said can be constructed with rj, 5 and A 

rrahcd ac / , icdcb . bd\ . i ■/ cb ^afcd cd Acfab\ 

nxy (adxOy + bdxdy + crjxyf] ) ~\- d i{r)fydxA -rjf^dyA-' ) (61) 

The last condition forces, that all four parameters a, b, c, d must be real. Every specific set of parameters (a, b, c, d) 
describes a different theory. Because the Lagrangian is a simple sum, it can also be written aj^ 

jC{a, b, c, d) = aCa + bCf, + cLc + ddd- (32) 

All individual terms Ca, ■ ■ ■ , are real for arbitrary complex r^^. Some of these Lagrangian terms are listed 
explicitly in the appendix lA. 21 Comparing (|32p with eq. (|28p . we see that the matrix Lagrangian is represented 
as Cz — C{1, —1, 0, —\)- In section [3731 is shown, that also the Lagrangian of the EiNSTElN-theory Ce — R can 
be expressed by this formula as £e ~ ^i^j 

Similar decompositions of the Lagrangian into a sum of terms, mostly in the teleparallel context, can 
be found in [36] and [26] (eq. (17) there). Also Itin [18], following the coframe description, gives a 3-term 
decomposion (eqs. 3.3 - 3.16) as most general form, which is for real r^^ equivalent to the first three terms in 
our eqs. H3ip - (|33p . Of course, none of them has a parity violating (PV) term ~ jC^, because for real tetrads 
obviously holds Cgj = 0. 

However, a similar term Cpy — rat"' with real tetrads was discussed in [29] as a possible cure for the initial value 
problem mentioned in sec 14. II But later it was shown that this term has to be rejected, because it leads to a ghost 
for the linearized theory ([21], P- 1219 and [30], p. 751). For the complex theory, presented here, the situation is 
quite different, because of the factor i the terms decouple (for all test cases with real tetrads), as demonstrated in 
sections 14.21 and 14.31 A deeper analysis of this, in connection with the discussion of the possibility of real tetrads, 
should be left to future work. 

Now we derive the geometric stress-energy tensor from the general form in eq. pO|l by variation of the tetrads 
(fe^. As usual, we consider and (e°)* as independent functions. Consequently we only have to variate rjj^^ and 
then express Sr^^^ in terms of So one gets as variation simply 

5C = {ryj*Hfy'=' = 5rl,uf = 1^,0^'^ . (33) 



We have defined here a new fundamental symbol f/"'' = {rl^)*HSy"'. It is a linear form of the ('"^^)* with 
constant coefficients. Because r^j, is antisymmetric with respect to the lower indices, only the antisymmetric 

part J/i"^' Ux^ — Ux°', which has also 24 components, is needed in (|33p. In the following is shown, how the 
stress-energy tensor is to compute using this symbol. Its explicit form, i.e. the form of the factors H^y'^'^ and the 
constants (a, 6, c, d) axe not needed for those general derivations. 



This is not the most general form for complex r^^. There exists e.g. a second, parity violating term 
~ '"aijC^cd)*^^!'^'''"^''' "which is not used here. 

Please, do not mix indices and parameters. Variable indices can never occur as factors. 

The summands Ca,- -,Cd are defined by eqs. II3UI I and Hijll l. and we note the correspondence Cd = ~2iCi from 
comparing it with the definition of d in eq. II28II . 

The variation of Se^ does not affect (r^f^)* because it is constructed (eJJ)* and their inverses only. On the other hand, 
the variation of the complex conjugated i5{e")* gives the same eqs. 
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3.1 Stress-energy tensor for the generalized Lagrangian 

For the computation of T from the Lagrangian C, the variation of all terms of the action integral must be 
expressed by the variations of covariant tetrads (fe^, so we need the formulas for the inverse (contravariant) 

tetrads and the absolute value of the tetrad determinant ||e|| '='' -^|e|iej* = y/—g, which are derived from the 
orthogonality relations: 



&^ = — eae"&Q and (5|e| = |e|e^&!^. (34) 



Inserting this, one gets 



a / u J a , u a 'y\ c h , c a / u a u oc\\ 
h / a "7 a 7\i£-a' / 1^ oc u 



^li'<'hfel - rlbe-Z) + Se^.aie'tet - e^e? ) (35) 



and the total variation of the action integral becomes 



51= I d^xdy\\e\\C)= I d''x((^|e||/:+ ||ell5/:) = / d'^x\\e\\[Sel;{^elC AD + Seif^aBl""). (36) 



The here introduced new expressions AJ^ and S^" ^'"^ compute by inserting the eq. (|35|l into 

eq. (|33|l . which expresses them by u'^^^ 



K = i^lf^l - rLeJWP = r^hfejiui' - U'J) = r^/e^t/i^'l (37) 

and the second is obviously the antisymmetric expression = —B'^: 

= (e)et - eJefWi' = eje^ui^'^ . (38) 

As usual, the variation term &^,q. in eq. (|36p is eliminated by partial integration (and neglecting the remaining 
surface integral) and this leads to the definition of the gravitational stress-energy tensor, here written as T^: 



51 = 1 d'^\\e\\5e'J,[^elC + Al- j^^iMBr),^] J d'x\\e\\ 5e'J; 



(39) 



with 

T;!.-lelC + Al^±-m\Br),c. . (40) 



This form with mixed- type indices (spacetime/tetrad, upper/lower) naturally arises from tetrad variation. If we 
want to transform it into a homogenous representation, we have to use a convention about the order of indices. 

Here we define T'^T =^ e'^'T^, i.e. the tetrad index should become the first. 

The above derivation of T from C is similar to the EiNSTElN-theory (Hilbert 1915), except that we 
used a more general Lagrangian and tetrads instead of the metric tensor. For the EiNSTEiN-case with 
jCe = ^(S/ii') 5/ii/,A) and 5jC = 5gfj,uTl^^^ it is easy to show, that one would obtain by tetrad variation like above, 
the tensor = ^hi^T^j^y which is equivalent. 

In the next section it will be shown, that a conservation law can be derived for the general stress-energy tensor 
defined in eq. (|40[) . that expresses energy-momentum conservation. 

However, for the general theory, in contrast to EiNSTElN-theory, where T^^^ ~ R'''^ — ^O^^R holds, the 

symmetry and reality of T^"^ is not guaranteed in all cases. This will be discussed in section |4l 
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3.2 Energy-momentum conservation 

In this section we derive a conservation law for the stress-energy tensor defined in eq. (|40p . As explained above, 
this definition holds for all gravitation theories, which are derived from a Lagrangian of the form H30[) . including 
Einstein- and matrix-theory. 

The easiest way to compute the covariant derivative is to use a tensor density (see e.g. 9 ), which here is 
defined by (the tetrad index "h" has to be transformed into a spacetime index "cr") 



77 =^ \\e\\e'iT;i = \\e\\{^S2C + eUl) - e'^iMB^'^),^. (41) 
We have to compute the divergence of this tensor densitvF^ 

77,^ = k\HC).a + i\\e\\eUl),^ - ei^,{\\e\\Br),^ - (Ilel|i3r),"7 (42) 



=0 

,Q7s h /|( 11 rTOn 



= ^{MC)^a + (lieiief,,„]B^^),7 - e^,^(||el|B,-["),o (43) 
= li\\e\\C).a + (||e||e|:,^„]B^^ - ei^MB^),! + 4nMe\\Br) (44) 



= ^(llell^).- - (||e||el,B^^),, = ^{\\e\\C),„ - c^,.7||e|lB,7 - e^,, (MB^''),^ (45) 

= ||e||(ieJ£+A°-T°) 

= llelK^^,. + liele'^^a + (e^??,.)*)/; - e^.^S"^ - el.(ie^£ + Af, - T^)) (46) 

= M\{\C,. + \{{ele'i<,ar - ele\„)C - et^^B^'' - 4,.(A^ - T^)) . (47) 

If we use the condition, that /I(e!^, e!j* , -J,, eQ*-y) does not explicitly depend on x^, we can compute its partial 
derivation with the definitions of the terms A, B 

'^''^ ^ ^e'^^'^'" ^ (^0,7),°- + cc- = + (^Z)*(e^,cr)* + -B^'^'fia.^o- + [B'^'^)* (e\^^„)* . (48) 

Inserting this in eq. H47|) . we compute the real part of the expression, where only one term on the rhs. remains 

^{J7^^) ^ WemeXaT^) ^ ^{eXaelT^) . (49) 

At last, we can easily show from the definition of the Christoffel symbols (considering only real tetrads), that 
for any symmetric tensor T^" = T"^ holds 

gXt^r^aT^" = eXc^hxT^", (50) 
so finally, if T is symmetric and real, the covariant derivative vanishes 

As conclusion it is to state, that the divergence of the stress-energy tensor is zero, if T is symmetric and real. 
This holds for all theories described by the Lagrangian of eq l|30|l . since the explicit structure of the symbol 
ul/^^ is not used in the above computation. 

We have to add, however, that for spaces which represent real matter distributions, the actual symmetry follows 
from the fact, that it equals the stress-energy tensor of matter Tj^J^ = T^^y This equation is usually derived by 
simply adding both Lagrangians and it postulates that matter acts as the source of the gravitational spacetime 
curvature. Esp. for the cosmological most relevant cases, the ideal fiuid approximation for matter is used, which 
is given by the real, symmetrical tensor 

Ti;:'^)^{p + P)u^u (52) 
where p is mass-energy density, p pressure and the 4- velocity. This is discussed in detail in section (4) 



Consider the antisymmetry of B'^ = —B'^ and the relation ejA^ = — eJeJ|,rjJj^_B"^ = (derived from eqn. 

0711 . II38I1 V Also used is the derivation of the tetrad determinant: |e|_CT = |e|e^e!^^ and in eq. II45II the definition of is 
inserted again. 

consider C = C* 
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3.3 Curvature scalar R and Einstein-Lagrangian in tetrad-form 

In this section we will show, that the Lagrangian of the EiNSTElN-theory can be written as special case of 
C{a, b, c, d) in eq. (|30|l . To prove this, we have to express the curvature scalar R by tetrads (we consider only real 
tetrads here). 

This is a quite lengthy computation, because one has to start with the complete RiEMANN-tensor, expressed by 
tetrads, and then to reduce it with R = r]"^^ Rmns- Similar computations can be found, with different notations, 
in various papers, e.g. [27] • Therefore, we have put it into the appendix IA.3I and will give here the result (again 

fb '''ba contracted form) 

T~> o I 1 c sb/o a , an x \ , xb /'e-oa 

7? = -2e rb,x + -rabV V-r^s + Vxcil rns) + n rxn . (53) 

This expression contains also second derivatives of the tetrads, namely the first term x. For the Lagrangian 
it is eliminated by partial integration: 

|e|e^Vb^A = - j {\e\e^^)^xn = - J v"'' iMes) ^xn = J \e\v^''rsrb , (54) 
so we finally get an Einstein-Hilbert-Lagrangian, which is bilinear in the first derivatives r^^: 

r sb 1 c sb/r, a . an x \ a c fQ/^bcd f cd 1 bd\ /cc\ 

Ce = V rsn - -r^i,r) {2rcs + Vxcil r^s) = rf^VgdV (SaSc ~ -^^^da - -VacV )■ (55) 

The same expression, in different notation, can be found in 24!, eq. (1) and [2S], eq. (6). Comparing this with 
eqs. dMll - dSll) gives 

CE = Cil,~^,-],0). (56) 



3.4 Isotropic coordinates and "viable" tetrad theories 

The term "viable" gravity theories is widely used in the literature. Nester [32] (introduction), defines it as "one- 
parameter class of teleparallel theories which agree with Einstein's theory to post-Newtonian order". Muench et. 
all [3T], give a similar definition of viable Lagrangians (p. 15), based on a three-parameter-set (ai,a2,a3), which 
is obviously equivalent to our set (a, b, c)0 

In this section we give a classification of tetrad theories defined by eq. ([32]). We show, that for all spacetimes, 
where isotropic coordinates can be used, a certain subset, described by the relation a -I- 6 -I- 2c = (including 
Einstein- and matrix-theory), have the same stress-energy tensor. Consequently, they have the same vacuum 
solutions, e.g. the fundamental SCHWARZSCHlLD-metric for spherical symmetry. Only those are considered as 
"viable" theories in the following sections. All others fail in the reality test. 

We request, that all viable theories must have real tetrads as solutions representing the SCHWARZSCHiLD-metric. 
Thus we can neglect the term dCj^, which is zero for real tetrads, in this section. (Its variation produces additional 
imaginary terms ~ i, however, which have to vanish independently, see section [42]) 

For concrete computations with tetrads, one must be careful not to mix the different index-types. Therefore we 
introduce here the symbol = as replacement term for the inverse tetrads, in this and the next sections. 

A static, isotropic coordinate system is defined with two real functions f{xi,..,x^),g{xi,..,X3) and the diagonal 
tetrads 

{eD = diag[f,g,g,g], {z^) = diag[-, -, -], \e\=fg^, (57) 

and it leads to the diagonal metric (g^jy) = diag[f^ , —g'^ , —g'^ , —g^], which includes the SCHWARZSCHlLD-metric. 
The not vanishing derivatives of the tetrads are 

eo,fc = /,fe, (^k,m=9.m, fc,m = l,2,3 

We now substitute / = exp{fi) and g = exp(A). The non- vanishing antisymmetric forms r^^ are (fc — fixed, all 
"r-triplets" - with three different indices - are zero): 

fe l^,k k k km ^.m , . 

'^Ofc = e[o,fc]^O^fc = — , rkm ^ ^[k,ni]ZkZm ^ —, k 771 



They give as viable class ai = l,a2 = —2, as = arbitrary (as = — i, for the teleparallel equivalent of Einstein's 
theory) . 
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To compute the stress-energy tensor, we need the terms Ul defined in eq. (|33[) . which are 

and the non-zero antisymmetric forms are (m = fixed) 

i[(a + 6 + 2c)/i^fe +2aA^fc] (58) 

-i[a^ -I- (2a + 6-1- 2c)A fe], m^k. 
9 

From the combination of the a&c-factors above, one can see, that all real tetrad theories with a-\-b + 2c = Q have 
the same [/-terms (up to a constant factor a, which we can set to a = 1, without loss of generality). 

Since the constants (a, b, c) appear nowhere else in the Lagrangian, those theories have the same stress- 
energy tensor. In section lOl it is shown, that the Einstein-Lagrangian is Ce = -C(l, — ^, — |, 0), which fulfills 
this criterion. Also matrix-theory Cz = C{1, —1,0, — ^) belongs to this class. Because a -I- 6 -I- 2c is the weight 
of all "r-doublet" -quadrats in eq. (|32p (terminology introduced in section eq. (llSp). this class is characterized 
by Lagrangians, which do not contain quadrats of r-doublets It is generated by setting 6 = —1 — 2c, which 
defines the set of "viable" theories by two real constants c, d 

Cv{c, d) -1 - 2c, c,d)= Ca- Cb + c [Cc - 2Cb) +dCd ■ (59) 

This class is investigated in the following sectional The value of the parameter c then defines the theory: c = (and 
d = — ^) describes matrix-theory and c = (and d = 0) is the Einstein theoryl3 For a convenient checking 
of the results, all terms of (|59|) and the [/-terms for this Lagrangian are explicitely listed in the appendix I A. 21 

In generalization of the eqs. H58|) . it can be shown (by the structure of Cx), that for tetrad fields, where all 
"r-triplets" r^^ =0, (x ^ y ^ z) are zero, the [/-terms of all viable theories are equal (independent of "c") and 
consequently the stress-energy tensor is equal to the EiNSTElN-tensor (except terms from of course). 

As the computations in the sections 14.3.11 - [4.3. 21 show, these viable theories also agree with the Einstein-theory 
in first and second PPN-order. 



=arfe-(2c + 6)r{;fe = - 
U]n^^ = arfe - (2c + 6)r"fe = - 



4 Comparison between Einstein- and matrix-theory 

For the comparison of different theories we use the general "viable" Lagrangian Ilv(c), defined in eq. H59|l above. 
For this we derive from eqs. pi|l and H33|) the following [/-terms, which are explicitly listed in the appendix eq. 
I|113p for convenient checking 

Ux — (V "x V "xPc + -^cji^q [Tcxj V \rcx) ) + .^cq 77 r]xy\r^^) 

+i(5St''* - SU''*) + iA^f'^SxSrl . (60) 

The computation of the stress-energy tensor for all test cases is then done with the following steps. 
(A) We start with the 16 covariant tetrads eJJ which represent the problem and compute (B) the determinant 
|e| and (C) the 16 inverse tetrads Cq, defined by the orthogonality eq. (D) compute the 24 coefficients 

'^bc ef/3,^]efe?. (E) compute the 24 [/i"''' with above eq. ^ resp. (fTT3)) . (F) compute C = ^r^b^t''^ 
and the 16 Al of eq. ^ and the 24 B^" of eq. (G) Finally compute with eq. (gOj and optionally 

yMT _ ri"^^emT^. These components of the stress-energy tensor then contain the parameter "c" and are valid for 
the class £v{c). 

For comparing the theories, we then have to use c = for matrix theory and c = — ^ (and formally set i = 0) 
for the EiNSTElN-theory. 

The above described computations are straightforward, but quite lengthy and error-prone. Existing software 
packages are either not well designed for these problems, or not free. 

That is why, we have developed "Symbolic" [3S], a small Java-program for such symbolic formula manipulations 
and the test of given solutions. It is a script-driven formula interpreter, especially designed for tensor calculus in 
GR, and produces TeX- and PDF-output files. It can be found, together with various sample scripts (nearly all 
test cases of this paper in tetrad formulation, as well as the corresponding problems for the Einstein theory and 
their results as PDF-files). Also available on this server is a web interface for testing it. 

^"^ We note, that the matrix-LAGRANGlAN can also be characterized as the only one, that contains no quadrats of 
"r-triplets" neither. 

The parameter d is not explicitly implemented, because it suffices to omit all terms ~ i (or formally set i = 0) to use 
d = 0. For d 0, its actual value plays no roll for real, symmetric matter tensors (e.g. vacuum) as it is demonstrated in the 
various test cases in section [1] However, the matrix Lagr. forces d = 
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4.1 "Unphysical" tetrads 

In the literature this kind of tetrads are discussed since the 1980-ies and by some authors they are considered as 
"death warrant" for the teleparallel theory (tetrad gravity). The first author, who presented them was Kopczyiiski. 
He showed in [20], that for some metrices the field equations are insufficient to determine the tetrads (resp. torsion 
tensor) completely. Then followed several papers, which tried to circumvent the problem, but all of them suffering 
from other serious physical problems 30 , 21]|5] . Esp. Nester [52] gives a very good overview about this dilemma. The 
essential statement of his paper is, however, that those tetrads are non-geneneric and occur only for very special 
solutions. Later work, using Dirac's constraint algorithm showed, that generic initial values have deterministic 
evolution while certain special initial configurations allow some undetermined evolution possibly only within a 
limited spatial region 

Here we show, that typical "strange" tetrads are excluded in the matrix theory, due to the parity violating 
term Cii in eq. A deeper, general analysis has to be done yet. A prototype for this kind of tetrads (compare 

[32| . p. 1008) is given with one arbitrary function xi^'^)' 



' cosh X sinh x ^ 



<^ I 1 I , I ^ 

^ sinh X cosh x / 

and it produces a flat MiNKOWSKl-metric — rj^i^, if x is reall^ The problem, that arose within previous 
tetrad gravity theories was, that the vacuum field equations = are identically fulfilled for any function x(^'') 
and thus do not pose any restriction on it. This fact obviously contradicts the assumption, that the tetrads (resp. 
torsion) possess a physical meaning, because they cannot be derived from some initial conditions. 
In the following we show, that for the matrix theory - also for solutions with real tetrads - there are non-vanishing 
terms ^ 0, accruing from C^, and thus this problem here does not exist. 

We have only non- vanishing 2 r-terms, namel\F^ 

^■03 = -X,ocoshx, J-oa = -X.osinhx- (62) 

For the EiNSTElN-theory it is obvious (the curvature tensor is zero), that the vacuum field equations are identically 
fulfilled. Since here all "r-triplets" are zero, it is consequently already clear from the considerations of section [57il 
that for real tetrads only the variation of jC^ can contribute to the field equations. 

For explicit computing, we do not list the intermediate (y-terms here (10 of them 7^ 0). We only note that 
£ = 0, all 16 AJ^ = and 10 -Bj" are 7^ 0. Finally, four do not vanish, which are explicitly 

tI = r| = i(x:oo + iXflfKe''-''' - e^'-^) + ix:oX,o(e^"^' + e^'-^) - (x:o)'e^"^' and 

= -T? = ^(x:oo + (x:o)')(e^"^' +e^'"^) + ^X:oX,o(e^-'^' -e>^*->^)-z(x:o)'e^-'^*. (63) 

We recognize from eq. (|63[) : 1. The constant "c" does not appear in any component of T, and consequently the 
real part of the stress-energy tensor is independent of "c", i.e. equal for all viable theories. 

2. For real tetrads (x = X*) follows — T2 = and the other two components T2 ~ i, Ti ~ i are only present 
for matrix-theory and the vacuum eqs. T2 = t1 = pose restrictions on xi.^'^) only here. 

3. The unique solution of = T2 = Q for general complex x (within the matrix-theory) is easily derived as 
simple linear function x{^^) ~ ^^'^ + c with two constants: k — real, but complex c. The free parameters c, k are 
then clearly determined by boundary conditions. 

As bottom line we resume again, that the problem solved in this section was not the existence of a solution, but 
the exclusion of physically unreasonable solutions. Of course, our computations here are no ultimate proof, that 
such solutions do not exist, but a strong argument. 



I.e. X = X*- For the matrix theory, however, we generally consider x(^'') complex valued function, which gives 
goo = —933 = 5R(e'<^~'<^ ), (?o3 = 0, from the definition II21II . 
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def 

For the inverse tetrads we use again the symbols = to distinguish them from the eJJ. They are given as 
23 = coshx, ^3~^o~~ sinhx- From this we get e.g. rjjg = — ^(Zyzl — ZqZ^) = x,0 coshx- 
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4.2 SCHWARZSCHlLD-solution 

In this section we show, that the important SCHWARZSCHlLD-metric is also a vacuum solution of the matrix field 
equations. From the considerations in section IXH it is clear, that the real parts of the stress-energy tensor are equal 
for all viable theories, i.e. also for matrix theory. It remains to clarify, however, that the additional imaginary 
terms do not pose unsolvable constraints. 

The tetrads to use are the same as in eq. (|57p . It shows, that it suffices to use only real tetrads for simplicity, 
i.e. real functions f{r),g{r). The computations (again following all steps from (A) to (G) on page 1121 which we 
do not list here) finally gives the following components of \e\T'^ (we list here 5 representatives, the other 11 are 
similar) 

|elro° = 2(341 + 3,22 + <7,33) - [g^ + g% + 9%)/9 (64) 
|elri° = 2j(/,2<7,3 - /,3ff,2)// 
leiTo^ = -4i(/,25,3 - ./,35,2)/5 

\e\Ti = {fg'^i - fg% - f9%)l9^ + (2/,i5',i + /5,22 + fg,33)/g + /,22 + /,33 
\e\Ti = -/42 + (/,2S,l - fg.l2 + /,l5,2)/5 + 2/343,2/3^ 

Their inspection shows, that imaginary terms ~ i only occur for and . They are zero for all spherically 
symmetric functions f{r),g{r), which was required. All other terms are real, and - since independent of "c" - 
equal for all theories. Hence it is obvious, that the vacuum solution is the well-kown SCHWARZSCHlLD-field. For 
completeness, we sketch some basic steps here. By the substitution g = we get for Tq — the simple second 
order eq. 

2(A" + -A') + A'2 = 0. (65) 

r 

This is solved by A' = — j.(^i^2r/M} ^^'^ leads to the well-known expression with an arbitrary constant ci: 

gkk = -5' = -e:cp(2A) = -c,{l + ^f. (66) 
The other components give two similar equations and finally lead to 

(1 - 

goo = /' = exp{2fi) = CO ^- . (67) 

which is the known metric for isotropic coordinates [2H] 1 [33 • According to the metric definition (|21|) the signature 
— , — , — ] is a forced result of the matrix theory. This is in contrast to other tetrad- or the Einstein-theory, where 
the signature must be postulated as additional assumption (e.g. as boundary condition for r — >■ 00). Unlike other 
tetrad theories, the matrix theory also does not presuppose the Minkowski metric, when the fundamental matrix 
Lagrangian eq. (Pl)) is considered. 



4.3 PPN-test 

In this section we perform a comparison between Einstein- and matrix-theory, based on the well-known PPN- 
scheme. It is shown, that both theories give identical results up to standard parametrized post-Newton (PPN) 
approximation order [21], |42|[^ 

4-3.1 Linear PN approximation 

For solving the linear field equations of the PPN scheme, we use a tetrad ansatz with 2x3 non-diagonal 
- generally complex valued - terms v^. ± /ij, {k = 1,2,3) and 3 equal space-diagonal elements e*^ = gSj, with 
f,g^ 1+0(6^) and h}^,v^. ~ 0{e^). It produces the metric tensor, which is used for the linear PPN approximation. 
Latin letters j,k,.. = 1,2,3 denote space indices. The symbols used here are in accordance to those for the 
SCHWARZSCHILD metric in eq. H57|) and section l472l because this ansatz can be considered as its generalization. 

eo = / = l + M. el=Vk+hk, eo=Vk-hk, e\ = g = 1 + \ 

We use tlie flat spacetime metric with the signature r] = [1, — 1, — 1, — 1], which has the opposite sign as in most 
GR-textbooks, and also for the metric results the opposite sign. Our form naturally evolves from the matrix theory (see eq. 
ill). It is also the form mostly used in relativistic quantum mechanics. 
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I f vi + hi V2 + h2 V3 + h^' 
vi — hi g 

V2- h2 g 

\v3- hs g 



(68) 



and gives the linearized metric 



500 


= 3?((eg)*e[j 


- (eo)*eo 


) ^ +l + 2K(At) 


= +1 + hoo 


= +1 + 


511 


= m^Tei 


/ 1\* 1 
- (ei) ei 


2K(A) 


= -1 + hii 


= -l + 0(e^ 


501 


= K((e[])*e? 


- (eo)*e} 


) ^ 23?(fei) 


def 

= hoi 


= 0(6^),... 


512 


^0, .... 











In this approximation only the real parts of /ij. enter into the metric, and the vj^ do not contribute at all. For 
the computation of the stress-energy tensor in linear approximation, as inverse tetrads Ca are to use the simple 
diagonals Ca ~ Sa and for the determinant |e| ~ fg'^ ~ 1. It shows in the following, that it suffices again to 
consider only real tetrads, i.e. all A, fi, h^.,v^. — real. This ansatz is valid up to the requested approximation order 
and solves all complex equations. 

We receive the following r-terms as step (D) in the general scheme (six representatives listed here): 

''01 = -(^ll +"l),0 +M,1 ri2 ^ hi^2 ~ h2A + yi,2 - V2A (69) 

''01 = {Vl ~ ftl),l - A_o ro2 = ('^l - /ll),2 '■12 = A, 2 ''23 = 

In the linear case, there is no need for the auxiliary terms Aj^, B'^" , since the stress-energy tensor can be 
computed from the ui^'^^ directly as ^ ~^h"a ~ ~^ha^' Here we need the components T^''' ^ rj^^T'^ which 
have to be equal to the ideal fluid tensor of the matter T^*/. For this comparison we compute the symmetrized 
y(M7) ^ yM7 _|_ y7M ^^d antisymmetrized components T^^'^l = T^'^ - T^^, for which we list six representatives 
here 

1^(00) ^ 2(Aai + A,22 + A,33) + O(e^) (70) 
T^'"' = /ll,22 + /il,33 - ^12,12 - /13,13 - 2A,oi + i(t;3,2 - '^2,3),0 + 0{^) (71) 



1^^(10) _ 4- ^ — — 9X 4- i(n,„ o — „ 4- n(t:^ 

' r(") = 2/12,02 + 2/13,03 + 2A,oo - A,22 - A,33 - ^,22 - M,33 + 2i(w2,3 - 1^3,2), 1 + 0{e^) (72) 



2 

r^^^^ = -^1,02 - h2fil + A, 12 + M,12 + i(^'3,l - '"1,3),1 + «(l'2,3 " 'y3,2),2 + 0{e'') (73) 



2 

_y[10l ^ _|_ l)(i;2,i2 + ■03,13 - ?;i,22 - «1,33) + i(/l3,2 ^ ft2,3),0 + 2i(l'3,2 - V2,?,),0 + 0{e^) (74) 
iy[21] ^ (-4^ _|_ l)(i;2 01 - t;i,02) + i(/ll,13 + /i2,23 + ^13,00 + ^3,33) + '(A, 03 " P'.Os) (75) 
+i(''3,00 + ''3,11 + 1'3,22 - t'3,33 - 2«i,i3 - 2li2,23) + C(e^) 
The inspection of these terms shows the following. 

1. The Einstein theory is given by c = — ^ and formally i = 0, which results in - as it must be - all rt'"^! = 0. 
The symmetric terms are simplified with the usual 4 gauge conditions, which read here 

A,fc + M,fc = 0(e^) and 2/ifc,fc -f 3A,o = O(e^) 

The fact, that the field variables v^. do not enter the field eqs., is a consequence of not contributing to the metric 
in this approximation. This is a basic problem of the - a priory symmetric - TEGR theory: the number of field 
variables exceeds the number of field eqs., resulting in free fields [27j . 

Since it is known, that the Einstein-theory gives the correct results, there is no need to compute it here. We 
sketch here only the basic computation steps, needed in the following, after 28,42 . The matter tensor of a perfect 
fluid is in this approximation given as (with rest-mass density p, pressure p and velocity itji^ 

= -Svrp, T-'" = -Sirpuj, T^'' = -8Tv{pujUk + pS^'') (76) 

The known result is 

l^=^hQ0 = -U + O{e^), h, = ^hoj = lVj + ^Wj+Oie^), \ = ~^hii=U + Oie^),... (77) 



The geometric stress-energy tensor T"^ in our notation equals — SttX the matter tensor T^P , and Uj = u 
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with the auxihary fields: 



2. For the matrix theory we have to show, that all additional terms in eqs. H7ip ... (|75[) (compared to 
the Einstein-case above) are zero in the requested order, and consequently the metric is the same. The above 
computation has made no use of the fields uj,, which thus can be freely chosen (within the order limit v^. ~ 0(6"^)). 
We set it here as with a simple ansatz from one potential "«" 

Vk=v,k- (79) 

With this ansatz all additional summands in the symmetric eqs. (|71[) ... (|73|) vanish. Also the term i(h3.2 — h2,3),o ~ 
0(6^) in ((Til does not contribute in the required order T^'^ ~ 0{e^) and finally the last eq. (ffSjPl 

ir[211 = t{h3^00 + "3,00 + ^A,03 - ^",3) « *(iA,03 - ^",3) = 0{e^) (80) 
leads to an additional "gauge" condition for the potential v 

Av = Ia^o- (81) 

It is solved with the help of the "super-potential" x(x, t) '^=' — f p(t, x')|x — x'ld'^x' (defined in [32], P- 94) and 
gives V = —^X,0 and finally 

As bottom line of this section we can state, that all T^" and consequently also the metric for the matrix 
theory are identical to those of the Einstein theory. 



4.3.2 Second PPN-order 

To perform the second order calculations we have to determine qqq up to 0{e'^). For this we have to use the same 
ansatz (|68p for (e) with more accurate inverse and the gauge conditions of section 14.3.11 already implemented 



(e) 



/ / + ^2 + ",2 ^3 + «,3 \ / 1// -ftl - ".1 -^2 - ",2 -^3 - ",3 " 

V.l~hi 1/f 

V.2-h2 1/f 

\v.3-h3 1/f 



1// 


-hi ~ v^i 


~/l2 — ",2 


~h3~ 


hi - v^i 


f 








h2 — U,2 





/ 





^13 - V.3 








/ 



with / = e^. The inverse tetrads (z) = (e) ^ above are sufficiently accurate up to O(e^) and also |e| = 1/ f'^ +0{e^) 
. Again, all tetrads can be considered as real, also for this approximation order. We skip all intermediate steps 
(D)...(G) here and give only the requested = r]"^'' z^Tjj^ . 
We obtain accurately up to order Ofe^lF^ 

= + A + t\3 - 2(ai,ii + M,22 + At,33) (83) 

+2j^i^3(/ii,2 - /i2,l) + 2i/i,2(ft3,l - ftl,3) + 2i/i,i(/i2,3 - ^13,2) + 0{e^) 
= Ma + t\2 + A - 2(/i,ii + M,22 + M.33) + 0(6^). (84) 

Again r'^'^ does not depend on the value of the parameter "c", i.e. it is identical for all theories of this class. It 
is consequently clear without computation, that matrix theory is up to this order not distinguishable from the 
Einstein theory. 



23 Consider ha, 00: ""3,00 ~ C(e^) 

2^ The terms 2ifi^:j(hi^2 — '12,1) -f ■ ■ ■ ~ C>{e'') can also be neglected, because ^ is needed only up to 0{e^). 
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4.4 New vacuum solutions 



The aim of this section is, to present some new vacuum solutions, which the EiNSTElN-theory does not possess. 
In the light of section \4A\ where we showed, that the matrix theory includes additional constraints, these extra 
degrees of freedom are not quite obvious. It might be possible, that solutions of this type can help to solve the 
galaxy rotation problem without the obscure "dark matter" , [51l41| . 

We use the following simple, static tetrads, with all diagonal elements — 1 and three real functions 
Vj^(x^ , x'^ , x^), k = 1,2,3. Here exact vacuum solutions can be computed quite easily, because |e| = 1 holds, 
and the inverse tetrads are also simple: 



(e) = {el) 



^ 1 Vi V2 V3' 

10 
10 

Vo 1 



f 1 —VI —V2 

1 
1 

Vo 



(85) 



The resulting metric is 

goo = 1, gok = gik = Vi^k - Sik- (86) 

The only three non-vanishing rj^^ terms are from the definition eq. (|18|) 



''12 = t^l,2 - ^'2,l, ''13 = ^1,3-^^3,1' ''23 = "2,3 - '^3,2 

and the stress-energy tensor can be computed exactly, following the remaining steps (E),.., (G), with 5 
repesentatives listed here (the 11 others similar): 

To = 3c(t;i^2 - W2,l)^ + 3c(t;i,3 - 1^3,1)^ -I- 3c(u2,3 - "3,2)^ 

+2c{vi^22 + 1^1,33 - '^2,12 - '^3,13)wi + 2c(l;2,ll + "2,33 - 'i'1,12 - W3,23)l'2 

+2c(u3^11 -I- 1^3^22 - '^1,13 - '^2,23)«3 
Ti = (1 + 2c)(u2,l2 + '^3,13 - '"1,22 - Vl,33) + 

c[{vi^2 - ■"2.\f' + (1^1,3 - "3,l)^ - («2,3 - 'VZ.2)^\V\ 

+2c(ui^3 - v-i,i){v2.;i - «3, 2)1^2 + 2c(?;i^2 - W2,l)('^3,2 - «2,3)"3 

+''^l(«2,3 - VZ.2),\ + ''^2(W3,1 - + ''^3('^1,2 - V2.\),\ 

To = 2c(-Ui^22 - ^'1,33 + ^^2,12 + '^3,13) 

Ti = -c(vx,2 ~ 'V2.\f' - c(wi,3 - ""3,1)^ + c(t;2,3 - "3,2)^ + '('^3,2 " «2,3),1 
T\ = 2c(u2,3 - «3,2)('^3,1 - «1,3) + '('^3,2 - «2,3),2 
We discuss here the vacuum solutions — 0. For the case c 7^ they force immediately 

Wl,2 - «2,1 = ^^1,3 - ^^3,1 = «2,3 - ^^3,2 = (87) 

This is exactly the condition for the flat (Minkowskian) spacetime r^^ = 0, which is consequently the only vacuum 
solution for the EiNSTElN-theory with c = — 1/4. 

For the matrix-theory we have instead c = 0, where the situation is quite different. Then the solutions of 

= are given by 

(«fc,j - Wj,fe),m = 0, i.e. Vk.j - Vj,k = Ckj = const., (fe, j, m = 1, 2, 3), (88) 

which is obviously a generalization of eq. (|87|) . The three antisymmetric constants c^m ~ ~emkj which define 
an axial vector, offer new deg rees of freedom as linear functions Vf^ — \^km-^ 

™, that in the EiNSTElN-theory all 

must be zero. We have to discuss however, that these solutions - although exact - cannot be regarded as global 
solutions, because the associated metric is not asymptotically flat. But they could probably be considered as 
regional approximation of similar generalized tetrads, which have to be found yet. 

In the realm, where the above solutions are approximately valid, they significantly modify e.g. the motion of test 
particles. This should be shortly sketched by a simple example. The relativistic equations of motion for small 
velocities vP ~ l,u^ <^ 1 give approximateljj^ 

^ (v^^k - Vk^m)u^ (89) 

If the motion then is considered as rotation inside a plane, perpendicular to the axis cj,„, we find a constant 
angular velocity, i.e. the tangential velocity is proportional to the distance from the axis. If we consider this as 
galaxy rotation, this increase is too fast, compared with the known flat rotation curves, [5ll41|. but this could be 
surely attributed to the simplicity of the tetrad ansatz eq. H85|) . 



The relevant CHRlSTOFFEL-symbols are \{gom,k - 9ofe,m) = \ {vm,k - Vk.m) 
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5 U{1) Noether-current 

Noether's theorem tells us, that every symmetry of the Lagrangian leads to a conserved current. The simplest 
case for matrix theory is the abelian (7(1) symmetry ej^ — >■ e^'^e'^, as explained on page|6l This current has no 
counterpart in real tetrad theories. In DiRAC's theory, however, it results in the conservation of charge [S]. 

To derive it here, we define the complex tensor /" '^^^ ||e]|B,'["e^ = \\e\Ku]f^^ (definition eq. (O). Inserting 
the field eqs. (gOj gives (with T T^eij)E!| 

1% = {\\e\\Bl°'e'^)^a = {\\e\\Bl°')^aey + ||el|B^"e!j_Q 

= Mi^elC + Al T^Je^, + l|ellB;[%> „ = M{-T + ^B^^^ef^^^^) 
= ||e||(-T + Z:). 

Since £ = real (and also T — real assumed) follows Ss(/").q — 0, and consequently is the conserved real (7(1)- 

Noether-current to define as the imaginary part J" = £s(/"). By inserting the (7-terms eq. (|60p for the matrix 
theory (with c = 0, see also (|28|l and (|113|l ) we get 

r = |lel|e?t//''I=||e||e?(-2^^=r*+3^^^*) 
and finally the conserved real current 

J =3(1 ) = -(/ -/ ) = llel||_2?7 (e^rc-Cb rc) + -{t e^ + e^ t )\. (90) 

Therein the first term vanishes for real tetrads (it contains only r-doublets) and the second term contains only 
r-triplets. 

The physical interpretation of this current is yet unclear. Its explicit computation shows, that it is zero for all 
astrophysical test cases in section|4l including PPN-tests up to order < O(e^). Hence it can clearly not be identified 
with a macroscopic matter flow. However, J° is not zero for the new vacuum solutions in section [4^ 



6 When are real tetrads possible? 

For comparison with existing "real tetrad theories" , we discuss here a modified matrix theory, which is described 
by the Lagrangian £(1,-1,0,0) = £a — Cf, in eq. ((32)) (i.e. without the PV-term Cij, resp. d = 0), and only 
considering real tetrads, for briefness labelled here as "real matrix theory" (RMT). Its i7-terms are represented 
by setting c = 0, and formally i = in eq. IjOOp . As explained in sec. 13.41 it belongs to the set of viable theories, 
which is widely discussed in the literature [311132] . The gravitational field eqs. of this "RMT" (e.g. for the vacuum 

= 0) are then a set of 16 real eqs. for the 16 real tetrad components e^. 
The "complex matrix theory" - presented in this paper - differs from this "RMT" by additional terms T^^^^ ~ * 

in the stress-energy tensoiPl in eq. (|40l) . which originate from the variation of £^ (as an example see the linear 
PPN-tensor in the eq-system ((TO]) ... (ffS)) '). 

If therein the tetrads are still constrained to be real, these terms are purely imaginary and decouple from the real 
parts T^j.-^^ (equal to the matter tensor) and build 16 homogeneous (non-linear, second order), partial differential 

eqs. T^.^l^ = 0, which are then additional and independent compared to the corresponding RMT. In this case, we 

consequently have a set of 32 independent real eqs0 for only 16 real tetrad components, which is expected to 
have generally no solution. 

However, one remarkable result of the test cases in section |4] was, that they all actually can be solved with real 
tetrads (for the PPN-tests they are real up to the required approximation order) 23 Therefore we shortly list 
the general form of these additional conditions in the following, although we are not yet able to give a complete 
mathematical and physical analysis of the solvability with real tetrads. 

26 consider e^ejj = 4 and A^eij = -2C 

For general complex tetrads these terms T^^j^ ~ * are not purely imaginary, nor are the T'^^j^ real. This is only the 
case for real tetrads. 

again look at the linear PPN example in sec 14.3.11 where actually all 32 eqs. are solved with real tetrads. 
2^ Obviously all statements in this section stay true, if we consider a constant, unitary transformation of all tetrads 
eJJ — > e'^eJJ, which also does not affect the metric. 
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If we consider only real tetrads, all terms rj^^, Tc, t"", |e|, . . . are also real, and the imaginary part of the symbol 
[7];"^' of eq. H60[l . which builds T'^^^^ becomes 

ul;^'l^^{5^,t'-5le)+^A^f^%,frc, i.e. [/{J = 3it^ and A = 0. (91) 
If we define the term T^^^^^ = ^7^-^(l)/i' following 16 conditions result, after some formula manipulations: 

miTi ■[/ a Acxfm , f m Acfab \ a Acftnh , a.,tn-\ ! ri /nn\ 

T(i)h = i[(rhfA, r)xa + -rab^^ rihf)rc ~ A ^ nhfr^a + ey^t ,a] = ^ (92) 

By contracting with it is possible to derive explicit formulas for = • ■ •. For the trace results the simple 
divergence-eq. 

T(i)m = + Cm*™) = |^ (! eje^t™) ,q = (93) 

For the important linear case ej^ « 5^ we receive the following approximation, which can be expressed by an 
antisymmetric "superpotential" FJl"^ = 

Fh = iVfheca+VfaecJi) with = -F/i,h = 0. (94) 

The trace vanishes identically T™^^_^ = 0, since t™j = 0. For the solvability of this system it is also important, 
that the identity -F/™fcm = ^ holds (because of the antisymmetry of F). 

It remains to clarify, for which matter tensors the eqs. (|92p resp. H94|) have no solutions, which also satisfy (|40|) . 
i.e. complex tetrads are actually required. 

7 Conclusions and outlook 

Here is presented a new classical theory of gravitation, which is in most test cases (SCHWARZSCHlLD-metric, post- 
Newtonian approximation), identical to the EiNSTElN-theory. But unlike other tetrad gravity theories, it does not 
exhibit some typical physical unreasonable vacuum solutions. 

It remains to clarify, if the correspondence of the symmetry groups of matrix-theory and standard electro- 
weak theory in particle physics SL{2) x U(l) D 5(7(2) x (7(1) is merely a pure coincidence, or if there are deeper 
connections between both. If the latter is the case, this would surely be worth of discussing in another paper. It 
should be possible to extend the global symmetry to a local one by introducing new gauge fields, likewise for the 
GSW-theory. But this is a quite complicated task, also needing a lot of new ideas. Also the issue of complex vs. 
real tetrads, requires further investigations. It should be clarified, in which cases real tetrad solutions are possible 
and how to interpret the possible imaginary parts physically. 

As shown in section f44l there exists a novel type of vacuum solutions, which are not present in the Einstein- 
theory. Although the sources of the field are not yet identified, these solutions have interesting properties regarding 
the galaxy rotation problem. To describe the sources of these solutions, it might be necessary to consider non- 
symmetrical stress-energy matter tensors. Einstein spent his last years searching for a non-symmetrical field 
theory [9], which was supposed to incorporate also electromagnetism, but without success. We know nowadays, 
however, that a classical field theory will not be able to answer all questions, because the wave- function in quantum 
mechanics cannot be regarded as physical field. 

Also the cosmological implications of the matrix-theory should be investigated. A very preliminary, first test 
with the simplest real tetrads, which produce the usual cosmological Robertson- Walker metrices, gives additional 
imaginary constraints, which force a spatially flat spacetime metric, i.e. k = 0. According to current astronomical 
knowledge, the matter density is nearly equal the critical density and does not allow the discrimination of k, so 
there is no contradiction. 

A remarkable, quite new perspective of the matrix-theory to spacetime geometry are the absolute matrices 
e.g. in These matrices are by definition invariant under all space-time transformations. 
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A Appendix 



A.l Matrix calculus 



Here wc want to list some formulas for matrix calculations, which are needed for the computations in sections 12. II and 12. 31 
Although quite elementary, they do not appear in most mathematical textbooks. 



a) For quadratic n X n matrices A, B, . . . of arbitrary dimension n holds the following. 
• Matrix factors inside the trace can be rotated cyclicallj[3 



r(ABc ■ ■ ■ X) = r(Bc ■ ■ ■ XA) = r(c ■ ■ ■ xab) = ■ ■ ■ . (95) 

• The trace of a hermitian matrix A^ = A is always a real number T{A) = real, and also of the product of two hermitian 
matrices T(AB) = real (using eq. I I95I I). But this generally does not hold for traces of more than two factors. 

• For the variation principle we need the following theorem: 

! ! 
The vanishing of the trace T(5kT) = for every variation matrix &c forces the matrix eq. T = 0. 

b) The rest of this section holds for 2x2 matrices only. 
We define for a matrix A = ("'^) a "bar" -operation ("adjunction") as the linear map A = 

• It is obviously interchangeable with hermitian adjugation (A)^ = (At), fulfills A = A and the evident equations with the 
identity matrix I: 

ab = bA, Aa = aA = |a|i, A + A = r(A)i, r(A) = r(A). (96) 

• The product of two matrices AB obeys no definite transformation rule under T-transformations defined in l lf 511 . but the 
"bar-alternating" product AB transforms in a definite manner as 

AB TATt T tBT = T(AB)r. (97) 



The same holds for products of more than 2 matrices. 

• As a special case of above, the trace of a bar-alternating matrix product with even number of factors is invariant under 
T-transformations, e.g. 

r(. .. ABCD...) =inD. (98) 

• If X, y,z, u are hermitian matrices, representing Minkowski spacetime vectors, the expressions 

F(x, y) ''^Z ^(xy - yx), V(x, y, z) ''^Z ^(xyz - zyx), V4(x, y, z, u) isT(xyzu) (99) 

are: F(x, y) = area (non-hermitian, traceless, 6 real comp.), V(x, y, z) = 3-volume (hermitian, 4 real comp.) and 
V4(x, y,z, u) = 4-volume (real scalar), respectively. All three expressions change the sign on odd permutations and vanish 
for linearly dependent vectors. 

c) Relations including the base-matrices t,i: 

• For every matrix A hold the three identities (to derive from the orthogonality and completeness of the basis) 



r(Af^)rf =r(Af'')T^ = 2A, r^AT^ = -2A, f At^ = 2ir(A) = 2(A A). (100) 

• For any non-singular basis (|r| ^ 0) and any index-combination a,/9,7 holds 

r°'f^r'' - T'fl^T°' = -2ie°''^^^Tx and e^,3.,xr°'f'^T' = Girx (101) 
where e is the completely antisymmetric tensor, with the scalar components e0l23 _ ^ 

^0123 = — It], . . .. These formulas allow an explicit computation of the contravariant- from the covariant matrices and vice 
versa. 

• To compute traces of products of PAULl-matrices, like in the eq. II26II an "index shifting" technique can be used, which is 
shortly sketched here. It is based on the orthogonality relations eq. ((Sj, which can also be written as am<^i +(7iam = "iVml^- 
We get e.g. 

T{umm<^'' ■■■)= r{{2-n^i - (7,a™.)(T" ■ ■ ■) = 2r;™^,r(<7" • ■ ■) - Tiaicr^a'' ■■■) = ■■■. (102) 

Using this technique multiple times, in combination with the symmetry relations eqs. II95I I and II96II gives the requested 
formulas. One example with 4 PAULl-matrices is the identity 

ir((7°CTV'=CT'*) = (r)''S'"' - + V'"'v'"') - i^"'""^, (103) 

where zi"'"^'* is the completely antisymmetric symbol, with A^^^'^^ = 1. 



The simple proof starts with 7~(AB) = 7~(BA), which follows e.g. from the component representation 7"(AB) 
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A. 2 Some explicit Lagrangian terms expressed by the symbols r^^ 

The following explicit expressions are included, to allow readers to check some formulas in this paper. They are computed 
with the help of a small computer program for symbolic computations "Symbolic" 1381 (see page [T2ll . but can be easily 
verified by hand. For uniqueness, the antisymmetric r^^ are always selected by the index combination b < c. Then the 
contracted terms of eq. I I27I1 are explicitly given as 

ro = r^i+r^2+r^3, n = -rgi + + rfg, ra = -rga - + rfg, rg = -rgg - r^g - (104) 

i" = — + — rj2, = —''23 ^ '"03 + '"02> *^ = ''13 + '"03 ~ ''Ol ' = ~''l2 ~ ''02 + ''Ol (105) 

First we list some terms of the general Lagrangian in eq. II32I I. 

= V"^"l~mr^ = — (''oi ~ ^12 ~ ''l3)(''01 ~ ^12 ~ ''l3 ) ~ (''02 + '"12 ~ ''23)(''02 + ''l2 ~ ''23) 

-('-03 + + '■i3)(''03 + '■W + ^23) + (■'•01 + ■'"02 + '"03)(''01* + ^02 + 4^) i^d) 

de_f rnn a , b _ _ 0* _ 0* _ 0. , / , 1 n 2* , , , 2 x 1* , / , 1 n 3. 
l-b — V ^mbVna) — ^OVOl "02*^02 '^03"o3 + l.'^12 + ^^02 + i '^12 + '^01 ^^^02 + i'^13 + '^03J'^01 

+ (-'"?3 + '"0l)'"03 + i^a + '■03)'-02* + (-'•23 + '•02)'-03 + '"Ol''-01* + (-'-02 + '"0l)'"S'2* (107) 

+ (~'"03 + '"0l)''l3 ~ '■i2''l2 ~ '"w^'lS + ('"23 ~ '"l3)''l2 + (~'"23 ~ '"l2)'"l3 + ''m'^oI 

+ 1 "^03 + '^02^23 ' 12*^12 + I '^13 +'^12^23 "23*^23 + '^03'^03 ^13"l3 '^23^23 

de_f mn „cd a , b N. _ _y 0. _ r, 0* _ r, 0* , r, 0. , p 0. , r, 0* , r, 1 1. 
L.C — rj r]abV ^mcVnd) — ^'^Ol'^Ol ^'^02'^02 ^"03"03 + ^'^12"l2 + ^'^13"l3 + ^"23*^23 + ^'^Ol"01 

+2''02''02 + 2''o3»'o3 — 2ri{2»'J| — 2rj^3ri{J — 2r23T-23 + 2rQjrg* + 2rQ2r§| + 2rQ3rQ3 — 2rj2r'i2 

~2r-j3rj3 — 2r23r23 + 2rgjrQj + 2rQ2»'o2 + 2''o3'"o3 — 2ri2»'i2 — 2r^3rfj — 2r23r23 (108) 

Cx = Cc ~ 2Cij = —2'r]abt'^t'"' = 2(rj2 + ^'02 ~ '"oi)(''l2 + '"02 ~ ''01) + 2(rj3 + Vq^ — 'rQi){ri^^ + rgg — rg*) 

+2(r^^3 + rg3 - rg2)(ri + r^s - ^2) - 2(r^3 - r?3 + rf^^rl^ - r^^ + rf^) (109) 
The Einstein-Lagrangian Ce reads explicitly (for real tetrads): 

CE = Ca- ^Ct - jCc = 2rl>^rl2 + 2(4^ - rl^)rf^ - 2r°^rl^ + 2{r^2 + '"D'^M " ^^As (HO) 
-2(r-83 + '•i3)^i3 - ^('■'12)' - ^(^?3)' - ^('"23)' + 2r^irg2 + 2(r-^i + ^2)43 + r^^r]^ 
~2^'"''2)'^ +'"13^03 ~ 2 '■''"3)^ g'-^^s)^ ~ i^02 +^12)^01 ~ g'-^'Ol'^^ +^23^03 ~ 2^*^03)^ 
+'-23'-?3 + ^irf.f - (rl, + r?3)ro\ - ^(4,f - {r^o3 + r°23)42 " ^(^§2)' + (rf, - rl^^.^ + ^(rf^f 
The two terms of the matrix Lagrangian Cz = Cr + iCi in eq. II28II are 

C, = Ca-Ct = irl2 + r?3)r0t + {-rl, + r^-^yO'^ + (-rJa - rl.,)rp, + (rl^ - rgi)r?2* (HI) 
+ (^■03 ~ ''0l)''l3 + (^03 ~ ^02)^23 + (""02 + ^03)^01 + ('"12 ~ '"0l)'"02 
+ (''13 ~ ^0l)'"03 + (~''"02 + ^m)*"!! + (~^03 ~ ^23)'"l3 + (^13 ~ ''l2)''23 
+ (~'"l2 ~ ^02)'"oi + ('"01 + ''03)''02 + (^23 ~ ''02)^03 + ('"Ol ~ '"13)^12 
+ ('"23 + ^12)'"l3 + (~'"03 ~ '"l3)''23 + (~''l3 ~ '"03)'"0* + ("'"23 ~ '"03)'"02 
+ (^01 + ^02)'-M + (-'"23 + ^13)^2 + (roi - rl2)rf3 + ('"02 + '"l'2)'"23 
= (^23 + ^03 — '"02)('"01 — ''12 " ''is) + (—''13 " ''03 + ''0l)(''02 + ''l2 " ''23) (112) 
+ (-'"23 + 'i3 - ''?2)('"01* + ''02* +4-3) + (^12 + '"02 " ''oi)(''o5 + ''« + '"!) 
+ ( — ''13 ~ ''23 ~ '"03)(''l2 + '''02 ~ '"01) + (''12 — ''23 + '"02)(''l3 + ''03 ~ '"01) 
+ (''12 + ''13 ~ '"01) (''23 + ''03 ~ '"02) + (''02 + ''03 + '"01) (''23 ~ ''13 + ''12) 

Ca consists solely of "r-doublets" , Cx solely of "r-triplets" . None of the Lagrangians Cr, Ci, Ce contains quadrats of r- 
doublets. In the terms Ca, Ci,, Cc, Cr, Ce r-doublets and r-triplcts do not mix, while d consists solely of mixed products. 
Cz = Cr + iCi does not contain quadrats of r-triplets. All Ca, C\,, Cc have even parity and only d has odd parity. 
The generalized "viable" LAGRANGIAN of section [3.41 has the form £^(c) = Cz -\- c{Cc — 2Cb) = Cz + cCx (definition in 
eq. II59II . but only for the special case d = —1/2, see footnote II8II . The antisymmetrized [/-terms, defined in eq. II33II and 
II6OII . for this Cv{c) are explictly (we list here 6 representatives, the other 18 symbols are similar) 

t^0°" = -^''23 - ^''03 + irl*2 + r\l + rll (113) 
= -(2c -f \)tI\ + (2c + l)rl2* + 2cr«2* + irH + ir^^ + ir^ 

^1 ' = *'"23 ~ *'"l3 + *''l2 + ''02 + ''03 

[/[021 = _(2c -f lyll + (2c + l)r02* + 2crj2* + i^l + ir\l + irl'^ 

[7^2] ^ .^0* + .^,1* _ .^3. _ ^0, ^ ^,3* 

[/[231 = _(2c + l)r?2* + (2c + l)rll - 2crll - irH - irH - ir^^, 
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A. 3 Computation of Riemann-, RiCCl-tensors and R with tetrads 

The aim of this section is to compute the R scalar with the tetrad- formahsm of section [2. II to enable its comparison with 
the Lagrangian of the matrix-theory, as presented in section [3] 
The RiEMANN-tensor is defined as the [Ai'J-antisymmetric expression 

T50" ■> p(7 I r^o" pa per per pa q(t qcj qit 

" V ' " V ' 

The CHRlSTOFFEL-symbols therein can be expressed by the tetrads (we only consider real tetrads here, because they suffice 
to describe RiEMANN-spacetime) using the standard formula 

We introduce the new symbols 7^^^ by transforming the upper index into tetrad type a —> s 

and with them the covariant tetrad derivative is defined as the expression (in contrast to the /^''-symbols, the are 
obviously tensors): 

G" e" = e" - r" e" = e" - F" fllfi) 

= ^(e[M..]+'='"(^-A"=[..a]+e-^[M.a]))- 

In the following we also will need their tetrad components, which are with the definitions in eq. II18II 

In some references these termes, which are by definition scalars, are titled "Ricci's coefficients of rotation". Then we can 
compute the second summand 

s%ux = r^.x + n^r^u (119) 
= {<n^),x + = e-r;,_;, + r'^A<,x + <Ca) (120) 

= <n^,X + n^{-^"«X + e?<^aA) (121) 



J> p6 



and we eelF^ 



<^;.,A - e"4n>^i<.x - Cx) = <irU,x - (122) 

ij.u\ — ^^nXu ~ ni'X (123) 




With this we can compute the tetrad components as 

mnl — <^me„e; n — e,„e„e; \y^^ )^ - y^^^XM + M^^aA ~ iixyau) 

= ^rn^n^t^^xy.X^fi^^' + ^xyi^^^i) ,X " G%y — G^yie^e^) ,v) + r^n^al ~ ^ml^an 

~^iymn,X ^nyml,v ^xy'^mf^n^l yy^i.L'^v ) ,X \ ) ,v ) ~'~ al ^miyan 

= '^lSmn,X ~ '^nGml,X + G%y<^m<^n<^l {^^ti. i'^Z,X ~ +'^i''^M,A ~ '^x'^fi.v) + ^mnGal ~ ^ml^^n 

e" 

= 4(^f5L„ - 5JiC;),A + G^xyis^n^rl, + sie-^i - 5fe^j - gj„c, + g=,r;^„ 

= e^(5f - Sf.G:^i),x + GxyiS^n-rli + Sle-^, - 5f - S^r^^ + 5f J 
= ^iSfGLr. - S^n5:^i),x + GxyiS^n^rli + Slie^^i - F^^) - Sf (e^^^ - J), 



The derivatives of contravariant tetrads are obtained from the orthogonality relations as ^ = — ej" eJJ'e' ^ 
^(M,.),A - ^(M,A),. = (4,- + - (e^,A + e^^).- = - ej^^, = ef^_^,_^ = -ef^^^j ef^_;,,_^ 
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so we have finally the tetrad representation of the RiEMANN tensor: 



Remarkable in this representation is the fact, that it is completely expressed by the Q and thus the r-terms, which in turn 
can be expressed by the p-tensor-matrix. 

From this we get by contracting over first and fourth index the tetrad components of the RiCCi tensor as 



Rmn — R'mns (125) 



and finally the R scalaiP^ 



= 2e^''<,,;, +g5,(r,-"rL +^""^(^3 +'7"'(w;^+'?scCd))) "^''''^ 
= -2e^''r6,A + G;?^i('?""rL + »?'""(C« + r?"^scCJ) + v'^'r^n 

->[a:i,] (a;i/) (i^!/) [xy] 
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